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Abstract 

Spatial correlation characteristics of the anomalous density matrix x. in a slab of nuclear matter 
with the Paris and Argonne vx% forces are calculated. A detailed comparison with predictions of 
the effective Gogny force is made. It is found that the two realistic forces lead to very close results 
which are qualitatively similar to those for the Gogny force. At the same time, the magnitude 
of n for realistic forces is essentially smaller than the one for the Gogny force. The correlation 
characteristics are practically independent of the magnitude of x and turn out to be quite close for 
the three kinds of the force. In particular, all of them predict a small value of the local correlation 
length at the surface of the slab and a big one, inside. These results are in agreement with those 
obtained recently by Pillet at al. for finite nuclei with the Gogny force. 
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I. INTRODUCTION 



The problem of surface nature of nuclear pairing has a long history, see the review paper 
[l| and Refs. therein. First, it was formulated in terms of the effective pairing interaction 
(EPI) entering the gap equation in a model space in which the pairing problem is usually 
considered. Within the self-consistent Finite Fermi System theory, the use of a natural 
density dependent ansatz for the EPI has resulted in a strong attraction at the nuclear 
surface, being rather small inside nuclei [2j. Similar conclusion was obtained in the ab initio 
calculation of the EPI [3J]. Later, the surface enhancement of the gap function A(X) was 
found by solving the gap equation for the complete Hilbert space in semi-infinite nuclear 
matter with the realistic Paris force in [4| and with the effective Gogny force in Similar 
conclusions were obtained for a nuclear slab in Ref. (d| where the gap equation was solved for 
both the types of NN force simultaneously. It was found that, although there is a quntitative 
difference between predictions of the two calculations, both of them show a pronounced 
maximum of A(X) at the surface of the slab, X = L, where 2L is the slab width, the effect 
being stronger for smaller values of L. In more detail, the gap equation for the nuclear slab 
] for the Paris force and in |8| for the Argonne v 18 force. It turned out that 
predictions of these two absolutely different kinds of realistic NN-force for the gap function 
agree with each other within 10%, both yielding the ratio A(X ~ L)/A(X = 0) ~ 2. 

Recently, Pillet et al. [9( investigated directly spatial properties of the anomalous density 
x which determines the space distribution of Cooper pairs. Calculations were carried out 
within the HFB approach with employing the D1S Gogny interaction [10] for a set of Sn, 
Ni, and Ca isotopes. It was shown that Cooper pairs in nuclei preferentially are located 
with small size (2 — 3 fm) in the surface region. The relevance of this phenomenon to two- 
nucleon transfer reactions was discussed. It should be mentioned that earlier the correlation 
properties of pairin g fo r specific nuclei were studied by Catara et. al. [ll], Ferreira et al. 
12]. Bertsch et al. [3], Hagino et al. [y| and yet in several works cited in A similar 
investigation has also been performed for T = pairing in dilute nuclear matter [15I ]. 

In this paper we carry out an analogous study for a nuclear slab with realistic NN force 
(the Paris and Argonne vi§ potentials) and the Gogny force. Our goal is to compare pre- 
dictions for the correlation pairing characteristics of the Gogny force and of realistic forces 
in order to analyze, to what extent the effect found in [9j is general and independent on 
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the specific choice of NN force. It should be mentioned that, with small modifications, the 
nuclear slab configuration may be used to describe the so-called "lasagna" phase of the inner 
crust of neutron stars. 



II. MAIN DEFINITIONS 

To make the comparison easier, let us recall the main definitions introduced in In a 
inhomogeneous system, the anomalous density matrix is defined as follows: 

x(ri,r 2 ) = ^Wi(riK(r 2 ), (1) 

i 

where Wj(r), t>j(r) are the Bogolyubov functions. For a spherical nucleus, it is convenient 
to go to relative and center of mass coordinates, r = ri — r 2 and R = (ri + r 2 )/2. In 9] 
the anomalous density matrix was studied in such a way that the probability distribution 
of Cooper pairs, |x(R, r)| 2 , was averaged over the angle between the vectors R and r, 

x 2 0R,r) = i- f |x(R,r)| 2 ^. (2) 

In particular, the space distribution of the pairing tensor \>c(R,r)\ 2 was analyzed. The 
probability distribution of pairing correlations, 

P(R,r)=R 2 r 2 x 2 (R,r), (3) 

was calculated in 9]. To avoid misunderstanding, this quantity is not normalized to unity. 
The coordinate dependent local correlation length was defined as 

_ (Jr 2 ^(R;r) (Pr) 1 ' 2 

At last, the locally normalized pairing tensor was considered in the form 

W(R,r)= r % {R : r \. , (5) 
J n z (R, r) r 2 dr 

as it enters the definition of the correlation length. 

Just as in js}, we consider a nuclear slab embedded into the Saxon- Woods potential U{x) 
symmetrical with respect to the point x = with potential well depth Uq = —50 MeV and 
diffuseness parameter of d — 0.65 fm typical for finite nuclei. The chemical potential is 
taken equal to fi = — 8 MeV. To compare our calculations with those of we fixed the 



thickness parameter of the slab as L = 6 fm to mimic nuclei of the tin region. We use the 
notation r = (s, x), where s is the two-dimensional vector in the plane perpendicular to the 
x-axis. The system under consideration is homogeneous in the s-plane, therefore one has 
x(ri,r 2 ) — > x(R, r) — ► x(X;x, s), with the obvious notation. The definition is then 
rewritten as follows: 

_ j(x 2 + s 2 )\^X;x,s)\ 2 d"r 

* {X) ~ f\x(X;x,s)\*d?r • W 

As far as the correlation properties in the x-direction and in the s-plane are essentially 
different, it looks reasonable to consider them separately, 

t 2 {X)=£{X)+g{X), (7) 

with the obvious notation. 

In the slab geometry, the angular averaging procedure similar to that in Eq. (T5J) is as 
follows: 

x 2 (A» = ^- J \x(X;x,s)\ 2 ^ 8(r 2 - x 2 - s 2 ) d 2 s dx . (8) 

It gives the distribution of the pairing tensor for a fixed value of the 3-dimensional relative 
distance r. There is another possibility, just to integrate over s: 



X 



(X,x) = / |x(X;a;,s)| 2 d 2 s . (9) 



It yields the quantity x 2 (X,x) which gives the distribution of the pairing tensor in natural 
for slab variables. For brevity, we use the same notation for the integrated anomalous density 
matrix as for the initial one and the one in Eq. (jEJ). The arguments should help to avoid 
misleading. Note that x(X, x) and x(X, r) have different dimensions. 

III. CALCULATION RESULTS 

Methods of solving the gap equation and the Bogolyubov equations for a nuclear slab 
are described in p| for the separable representation of the Paris potential and in [8], for the 
Argonne v 18 force. The latter could be used for arbitrary NN-potential, and we repeated 
all the calculations of 8] for the Gogny force. To begin the comparison, let us start from 
infinite nuclear matter. The dependence of the gap on the density of nuclear matter for the 
three kinds of the NN-force is displayed in Fig. [TJ The correlation length (Tj0) in infinite 
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FIG. 1: The gap A(k = kp) in infinite nuclear matter 
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FIG. 2: The correlation pairing length in infinite nuclear matter 



matter can be easily found in the momentum space: 
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e = 



d_ 

dk 



x{k) 



d 3 k 



f\x(k)\ 2 d 3 k ' 



(10) 



Let us substitute in this equation x{k) = A(k)/2Ek, where Ef. = a/ (ifc — £ F ) 2 + A 2 (k), 
£fe = k 2 /2m*, Ef = k F /2m* . The functions inside the integrals both in the numerator and 
denominator of this relation are very peaked in the vicinity of k = k F and rapidly vanish 
outside the interval \k — kp\ < k F (A F /ep), A F = A(k F ). Usually one deals with the limit 
(Af/£f) << 1. In this case, one can substitute A(&;) = A F in Eq. (flQj) and evaluate the 
integrals analytically. A simple calculation yields 

where t> F = k F /m*. 

The correlation length for the three kinds of force under discussion found numerically 
from Eq. ( flOl) . with m* = m, are displayed in Fig. [3 For comparison, the approximate £ 
from Eq. (TTTT) with Argonne force is also displayed. It is seen that the approximate formula 
works sufficiently well in all the interval of fc F . Even at the maximum of the gap the deviation 
from the numerical result is of the order of 15%. 

One can see that at small density, k F < 0.5 fm _1 , the results for all three forces practically 
coincide. This is not strange. Indeed, although the Gogny force is an effective one, in the 
1 Sq channel under consideration it describes the free NN-scattering perfectly well for small 
energy values which are only important in this density interval. In this sense, the Gogny 
force could be considered as a semi-realistic force. Two realistic forces lead to close results 
for all density values. In the vicinity of the gap maximum, the difference between A Arg 
and A Par does not exceed 10%, and only at k F ~ 1.4 fm _1 , where the gap value itself 
becomes very small, the relative difference becomes larger. As to the Gogny force, at the 
density region kp ~ 1 fm _1 it leads to the gap values which are bigger by approximately 
(20 -T- 30)% than those for realistic forces. Correspondingly, the correlation length for the 
Gogny force is quite close to that of the realistic forces till — 0.8 fm _1 , and only at 
k F ~ 1.2 fm _1 the difference becomes large. The density dependence of the correlation 
length, £(/c F ), is qualitatively similar for all the three types of force. It consists of a plateau 
at 0.3 < k F < 1 fm _1 and two intervals of sharp growth, at k F < 0.3 fm _1 and k F > 1 fm _1 . 
In the latter, the value of £ Gog (&F) is growing with /cf much slower than that of ^ Arg (^F) 
and £ Par (/cp)- Note that at k-p ^ 1.2 fm -1 the difference between £ Arg (fc F ) an d £ Par (&F) also 
becomes noticeable. This is a manifestation of their behavior near the critical point k F at 
which the gap vanishes and transition to the normal phase of nuclear matter occurs. The 
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FIG. 3: The Fermi averaged gap Ap(X) 
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FIG. 4: The Fermi averaged effective pairing interaction V^ S (X) 

values of k F for the Argonne force and the Paris potential are a little different. This results 
in different behavior of £ Arg (/cp) and £ Par (/cp) in the region of kp ~ 1.5 fm" 1 . 

Let us now turn to the slab system. Before analyzing the correlation characteristics, it is 
instructive to briefly compare the EPI and the gap itself found with the realistic forces and 
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FIG. 5: The correlation pairing length £(X) in a slab of nuclear matter 
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FIG. 6: The correlation pairing length £ X (X) in the x-direction for a slab of nuclear matter 
the Gogny force. The "Fermi averaged" gap is displayed in Figj3l It is defined as follows: 

A F (X) = A(X,k F (X)), (12) 



where the local Fermi momentum is k F (X) = ^2m(ij,—U(X)) U(X)). This quantity 

characterizes the gap on average [l[ Q, Is]. We see, first, that all three functions Ap(Jf) 



S 



have pronounced maxima at X ~ L = 6 fm. The ratio Ap(X ~ L)/Af(X = 0) ~ 2 



for realistic forces and ~ 1.5 for the Gogny force, in agreement with jg]. Second, the gap 
Ap° g is significantly bigger than Ap ar and Ap rg , by approximately a factor one and a half 
at the surface and two inside the slab. It is worth to discuss this point in more detail. 
Let us consider 120 Sn as a "reference nucleus" . Its empirical gap value is estimated usually 

17]. Diagonal matrix elements of the gap found in ?], Q] 



as A ~ 1.3 MeV 
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for a slab 



with L = 6 fm are about 1 MeV which agrees with the above value, leaving ajroom about 
20-30% for the surface vibration contribution. The 
which is, in our opinion, too much, see discussion in 



atter was estimated in [17J as ~ 50% 
ifl. We consider the estimation of 18] 



at ~ 30% as more realistic, but, evidently, also too big, because of disregarding so-called 



tadpole diagrams 19j. Calculations of ]9] for this nucleus gave A ~ 2 MeV which, in our 
opinion, is too much, especially if to take into account that some additional contribution of 
surface vibrations to the mean field theory value of A should be! Thus, our observation in 
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ilations that the Gogny force overestimates the gap value agrees essentially with 

results of 

To understand the physical reason for the surface enhancement of the pairing gap with 
each NN force under consideration and bigger values of the gap for Gogny force, it is useful 
to calculate the EPI which we use in the two-step method of solving the gap equation [l|. 
Let us remind how this quantity is defined. In a symbolic form, the microscopic gap equation 
reads: 

A = VA S A, (13) 

where V is the free iViV-potential and A s = GG S stands for the two-particle propagator in 
the superfluid system. Here G and G s are the one-particle Green functions without and 
with pairing effects taken into account, respectively. In Eq. ( |T3l) . as usual, integration over 
intermediate coordinates and summation over spin variables is understood. Let us now split 
the complete Hilbert space S of two-particle states into two parts, S = S + S'. The first 
one is the model subspace So in which the gap equation is considered, and the other is the 
complementary subspace S'. They are separated by the energy E in such a way that S 
involves all the two-particle states (A, A') with the single-particle energies E\,Ey < E . The 
complementary subspace S' involves the two-particle states for which one of the energies 
£\,£\> or both of them are greater than Eq. Therefore, pairing effects can be neglected in 
S' if Eq is sufficiently large. The validity of inequality A 2 /(E — /i) 2 <C 1 is the criterium 



of such approximation. Correspondingly, the two-particle propagator is represented as the 
sum A s = Aq + A'. Here we already neglected the superfluid effects in the S"-subspace and 
omitted the superscript "s" in the second term. The gap equation (fTBj) can be rewritten in 
the model subspace, 

A = V? S A S A, (14) 
where the EPI should be found in the supplementary subspace, 

V^ = V + VA'V: S . (15) 

Note that the last equation has a strong similarity with the Bethe-Goldstone equation. 

As the analysis showed S|, the optimal choice of splitting corresponds to Eq = 
15 -7- 20 MeV. In a slab system, the EPI is calculated in the mixed coordinate-momentum 

r B.„,,...,,„ M „...,. m ..„„,,..^ 

form [3, LSI] with the Fermi averaged strength 

V c F ff (X) = J ^ F M F (I) ; I + ^-^). (16) 

The Fermi averaged EPI for the Argonne v 18 force and the Gogny D1S force calculated for 
Eq = 15 MeV are displayed in FigHJ We did not display the EPI for the Paris potential as 
it practically coincides with that for the Argonne force. One can see that both the curves 
behave in a similar way changing from quite week attraction inside the slab to very strong 
one outside. The reason for the latter is that in the asymptotic region X > L the V^ s value 
tends to the quantity which is very close to the free T-matrix taken at the negative energy 
E = 2/i. To be precise, the limit is equal to T(E = 2\i) if the separating energy Eq = 0. In 
the case of Eq ^ the limit is equal to some "T'-matrix" which is obtained by solving the 
same Lippman-Schwinger equation as the T-matrix, but in a cut momentum space, because 
the contribution of nucleons with total energy less than Eq must be pulled out. As it is 
known, the Gogny force leads to the scattering length in the 1 Sq channel a ~ 12 fm Q]- 
It differs, of course, from the experimental value of a ~ 18 fm which is reproduced by any 
realistic force, but not so much. In any case, the virtual pole of the T-matrix for the Gogny 
force is close to zero as it should be. Therefore the analytical continuation of the T-matrix 
(or T'-matrix) to rather small negative energy E — 2{i results in an enhancement of the T'- 
matrix (~ —950 MeV fm 3 ) in comparison with typical values. This enhancement is not so 
strong as in the case of the Argonne force (~ —1100 MeV fm 3 ), but it is equally significant. 
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FIG. 7: The pairing tensor distribution x 2 (X,x) 

The inner value of the EPI for the Gogny force is quite small (~ —160 MeV fm 3 ), but bigger 
;han for the Argonne one (~ —95 MeV fm 3 ). In our previous study with realistic forces 
l|, lD, |8] we explained the surface enhancement of the gap in terms of the sharp variation of 
the EPI at the surface. In the case of the Argonne force, the ratio V^/V"^ ~ 12. For the 
Gogny force, it is about 6. This is also a big number leading to a surface enhancement of 
the gap, but not so pronounced as for realistic force. FigfJ] explains also why the Gogny gap 
is so big. It is well known that the pairing gap depends on the interaction strength in an 
exponential way. In [3] it was found that 1% variation of V^j leads to 5% variation of the 
gap. It explains why the gap function for the Gogny force is in 1.5 4- 2 times greater than 
the one for realistic forces. 

Let us turn now to the correlation pairing characteristics. The total correlation lengths 
£ and the correlation lengths in x-direction £ x found from Eq. (J7J) for each kind of force 
are displayed in Fig. [5] and Fig. [6j correspondingly. The curve £ Arg (X) is quite similar to 
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£ Par (A), the difference is about 5%. Both have a minimum shifted a little from the surface, 
X = L, in direction to the free space, both have a pronounced maximum in the vicinity of 
X = and grow rapidly to the right from the minimum. Qualitatively, the curve £ Gog (X) 
behaves in a similar way, but the maximum value at X = is approximately two times 
less than those for realistic forces. The minimum of £ Gog (X) is also shifted from the point 
X = L, but the value of the shift is less. Such a behavior of each curve £(X) qualitatively 
agrees with naive local density approximation (LDA) predictions. Indeed, inside the slab, 
the local Fermi momentum kp(X ~ 0) ~ ^/2m(/x — Uq) — 1.4 fm" 1 which corresponds to 
big values of £ in infinite matter (see Fig. |2J). The same is true at X > 8 fm where the 
local Fermi momentum vanishes. But under more detailed examination, deviations from the 
LDA predictions are significant. To illustrate this point, we display in Fig. [5] with a thick 
line, in arbitrary units, the density distribution p(X). Within the LDA, the correlation 
length £(X) should show a plateau inside the slab, just as the density does. Instead, £(X) is 
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FIG. 9: The probability distribution r x (X,r 



decreasing rapidly with increase of X till X ~ L. Qualitatively, the coordinate dependence 
of the function £(X) in the slab reminds us that of £(R) in spherical nuclei found in 
Both of them have minima at the surface region and pronounced maxima in the center. But 
any quantitative comparison is hardly possible since, for the problem under consideration, 
the properties of the two systems are essentially different. Indeed, in a spherical nucleus all 
particles move in a finite space limited by the nuclear surface. On the contrary, in a slab the 
particle motion is limited only in the x-direction. In the s-plane, the motion is free which 
leads to very big values of £ s in Eq. © and £ close to that in infinite system. As to the 
£ X (X) function, it should be much closer to £(R) of 9). As it is seen in Fig. [BJ two curves 
^ Ig (X) and £^ ar (X) practically coincide. Deviation of £^ og (X) from the both is much less 
than in the case of the total correlation length. It doesn't exceed 15%. All three curves have 
common minimum at X ~ L, with the value of £™ in ~ 1.5 fm. It is not far from the value of 
£imn ~ 2 fm found in [9] for Sn isotopes. Evidently, the difference is mainly due to geometry 
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FIG. 10: The profile functions r 2 x 2 (X = X ,r) 



effects. Thus, for a nuclear slab, the correlation length of pairing in the x-direction at the 
surface, calculated with realistic and semi-realistic Gogny forces, is very small, in agreement 
with conclusions of [9|. 

To visualize the pairing tensor distribution, we draw a 3-dimensional plot in FigJT] for 
the >c 2 (X,x) function given by Eq. for the Argonne force. We see that there is a set 
of maxima at x — 0, the highest one being near to the slab surface, X ~ L. FigJS] shows 
the profile functions x 2 (X = Xo,x) for several values of Xo corresponding to the maximum 
positions. The nearest to the surface maximum there is at Xo ~ 5 fm, the neighboring 
one, at Xo ~ 3 fm. There is also a pronounced maximum at Xo = 0. The surface peak 
is very narrow, in correspondence with FigEJ On the contrary, in the case of X = a 
comparatively sharp peak is accompanied by a flat base plate which results in a big value 
of £ x . The similar profile functions are drawn, in the same figure, for the Paris and Gogny 
force. The Paris curves are again quite similar to the Argonne ones. As to the Gogny force, 
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FIG. 11: The locally normalized pairing tensor W(X,r) 

the similarity takes place only at qualitative level, absolute values of x 2 being bigger than 
those for realistic forces. This is the result of bigger values of the gap itself for the Gogny 
force, as illustrated in Fig. [31 We see that Ap ar is a little bigger than Ap rg , therefore (x Par ) 2 
is, on average, bigger than (x Ar «) 2 . On the other side, Ap° g is significantly bigger than the 
gap for realistic forces. As the result, (x Gog ) 2 significantly exceeds the microscopic values 
as well. Note that the correlation length £,(X), Eq. (jSJ), does not depend on the magnitude 
of x, giving rise to a much smaller difference between the Gogny and realistic forces. The 
same is true for £ X (X), Eq. (J7J). The surface peak dominates, to some extent, for all three 
kinds of forces, and this effect for realistic forces is even stronger than for Gogny force. 

The x 2 (X, r) function given by Eq. (jSJ) multiplied by r 2 gives the probability distribution 
of paring correlations similar to Eq. in the case of the spherical symmetry. It is displayed 
in Fig El again for the Argonne force. Now the main maximum positions are shifted from 
r = 0tor~l^-2fm due to the factor r 2 . The surface maximum at X ~ L is even 
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FIG. 12: The profile functions W{X ,r) 

more pronounced than in FigJTl To compare results obtained for different forces under 
consideration, we again draw the profile functions r 2 x 2 (X = X , r), see FigfTUl Again, just 
as in FigJHl the Gogny force results are significantly bigger than those of realistic ones in 
the magnitude but are very similar in the form. Again the surface maxima dominates and 
again the surface enhancement is stronger for realistic forces. 

To make the comparison with [9] more complete, we display in FigJTl] the locally normal- 
ized pairing tensor which is defined as 

r 2 x 2 (X, r) 



W{X,r) 



(17) 



J >c 2 {X, r ) d 3 r 

similar to the definition (j3J) in spherical systems. The profile functions W(X = X Q , r) are 
displayed in Fig{12] which is analogous to FigJH] and FigfTOj In this case, the Gogny force 
curves are close to those for realistic forces not only in the form but also in the magnitude. 
This is happens because the x 2 quantity comes in both the numerator and the denominator 
of Eq. (117p . the result being almost independent of the magnitude of h. Absolutely similar 
situation occurs with calculations of the correlation lengths £ and £ x . 
The denominator of Eq. (TTJ 



p(X) 



x 2 {X,r) d 3 r = / x 2 (X,x) dx 



(18) 



gives the total probability distribution of Cooper pairs integrated over relative coordinates. 
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FIG. 13: The integrated Cooper pair probability distribution p(X) 

Note that this quantity displayed in Fig{T3l just as P(R,r), Eq. Q, is not normalized to 
unity. Again we see that the result for the Gogny force behaves qualitatively similar to 
those for the realistic forces but its magnitude is significantly bigger. And for each force 
under consideration this quantity in the slab does not exhibit any surface enhancement. It 
occurred because, although the surface maximum for any force in FigJH] is higher than the 
central one, the correlation length in the x-direction at the surface is much smaller than 
inside the slab that makes the integral over x to be smaller. An analogous effect should 
take place in nuclei, too, but in this case the "geometrical" factor R 2 in Eq. ([31) would help 
to survive to the surface enhancement in the Cooper pairs distribution if the probability 
P(R,r), Eq. ([3]), is integrated over the relative coordinates. 



IV. CONCLUSION 



Spatial correlation properties of nuclear pairing in a nuclear slab are studied with two 
realistic NN potentials, the Paris and Argonne v 18 , and with the phenomenological Gogny 
D1S interaction. The results obtained with the two realistic forces agree with each other 
with an accuracy of about 10%. But they agree only qualitatively with those of the Gogny 
force. The gap value in the slab for the Gogny force exceeds that from the realistic forces 
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by a factor about 2 , which results in rather bigger values of the anomalous density matrix 
x. Nevertheless, some of main conclusions of [9( obtained for finite nuclei with the Gogny 
force are confirmed qualitatively, especially the dependence of the correlation length on the 
position of the c. m. of a Cooper pair. This quantity does not depend practically on the 
absolute value of x, only the space distribution of the pairing tensor x 2 (X, x) being impor- 
tant. At the surface of the slab the local value of the correlation length in the x-direction is 
very small, £ X (X ~ L) ~ 1.5 fm, for all three kinds of the force under consideration. Inside 
the slab £ x becomes very large, i.e. of the order of the slab width or even more. Thus, in 
this point our results completely confirm those of 

The pairing tensor x 2 (X, x = 0) has several maxima, among them the ones at X = 
and X ~ L are most pronounced. And the one at the surface is a bit higher, especially for 
realistic forces. In this sense, we can speak of surface enhancement of the Cooper pair distri- 
bution. However, the total probability p(X) for a pair to have the cm. coordinate X, which 
is obtained by integrating x 2 (X,x) over relative coordinate, has no surface enhancement. 
Thus, the second conclusion of |9|] that Cooper pairs in nuclei prefer to be concentrated in 
the vicinity of the surface should not be drawn for a slab. We explain this with different 
geometrical properties of two systems under comparison, with different "surface-to-volume 
ratio" in a sphere and in a slab. However, all surface enhancement features found for real- 
istic forces are qualitatively reproduced with the Gogny force. We trace this effect to the 
"semi-realistic" nature of the Gogny force which describes the low-energy NN-scattering in 
the 1 S'o-channel sufficiently well. It seems reasonable to suppose that all main conclusions of 
jsl will be confirmed qualitatively if the Gogny force in the pairing channel will be replaced 
by a realistic NN-potential. 
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